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We study Newtonian cosmological perturbation theory from a field theoretical point of view. We 
derive a path integral representation for the cosmological evolution of stochastic fluctuations. Our 
main result is the closed form of the generating functional valid for any initial statistics. Moreover, 
we extend the renormalization group method proposed by Mataresse and Pietroni to the case of 
primordial non-Gaussian density and velocity fluctuations. As an application, we calculate the 
nonlinear propagator and examine how the non-Gaussianity affects the memory of cosmic fields 
to their initial conditions. It turns out that the non-Gaussianity affect the nonlinear propagator. 
In the case of positive skewness, the onset of the nonlinearity is advanced with a given comoving 
wavenumber. On the other hand, the negative skewness gives the opposite result. 

PACS numbers: 98.80.Cq, 98.80.Hw 



I. INTRODUCTION 



The large scale structure in the universe has evolved from primordial fluctuations according to the gravitational 
instability. In the standard scenario of the structure formation, the primordial fluctuations are created quantum 
O |. mechanically during the inflationary stage in the early universe. After exiting the horizon, the fluctuations are 
evolved linearly; which is well described by relativistic linear perturbation theory. Eventually the fluctuations re-enter 
into the horizon. After that, it is sufficient to treat the evolution of fluctuations by means of Newtonian gravity. Due 
to the Jeans instability, at some point, the density fluctuations become nonlinear. In this stage, usually we resort to 
the N-body simulations. However, since the numerical simulations are time consuming, the analytical calculation of 
the nonlinear evolution is still desired. The standard perturbative expansion method is developed for this purpose. 
| In the quasi- nonlinear regime, the perturbative approach was successful [H 0, S 0, H S 0| • To obtain more accurate 
. results, however, the non-perturbative analytic method would be necessary. 
t^J- ' Recently, Crocce and Scoccimarro have developed a new formalism to study the large scale structure Q. They 
described the perturbative solution by Feynman diagrams and identified three fundamental objects: the initial con- 
^sO • ditions, the vertex, and the propagator. They have found that the renormalization of the propagator is the most 
important one. Based on this finding, they have observed that, due to the rapid fall off of the nonlinear propagator, 
' the memory of the cosmic fields to their initial conditions will be lost soon in the nonlinear regime. 

Following their work, Matarrese and Pietroni reformulated the cosmological perturbation theory from the path 
integral point of view and developed the renormalization group (RG) techniques in cosmology [Tl], 03 ( see @ for a 
slightly different approach). Matarrese and Pietroni have app lied their formalism to the baryon acoustic oscillations 
(BAO) which takes place around the scale k ~ O.lMpc -1 [l3j . On these scales, the nonlinear effects are relevant [id], 
Ha ]. They have found that the renormalization group method is useful to predict the BAO feature. Crocce and 
Scoccimarro have used their graphical approach to discuss BAO and found the renormalized perturbation approach 
gives a good agreements with results of numerical simulations [la ]. This result is further confirmed by Nishimichi et 
al. [13. 

These authors have discussed only the Gaussian initial conditions. Recently, however, it has been realized that the 
primordial non-Gaussianity can be produced in the inflationary scenario. If so, it is important to give a renormalization 
group formalism for the non-Gaussian initial conditions. Conventionary, the non-Gaussian curvature perturbation $ 
is characterized by the following form 

$(x) = $ s (x)+/^($2( x )_($2 ))j (1) 

where is the Gaussian field and /jvl is the parameter to represent a deviation from the Gaussianity. There are 
several possible observational tests to constrain the non- Gaussianity [H HE [H, HE HE HE HE HE Hj . The most 
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stringent limit comes from WMAP and the result is —58 < Jnl < 134 [2(J. Planck and other tests will give a more 
stringent limit or detect the non-Gaussianity [2l| . 

The purpose of this paper is therefore to extend the analysis by Matarrese and Pietroni to the non-Gaussian 
initial conditions. Starting from the generating functional of the multi-point functions, we derive the path integral 
representation of the cosmic evolution of the cosmic fields. In contrast to the previous work, the non-Gaussianity is 
incorporated into the field theoretical scheme. In particular, we obtain the formula for the generating functional which 
allows us to use the Feynman diagram method to calculate various statistical quantities characterizing the large scale 
structure in the universe. We also derive the RG equation for the effective action. As an application, we calculate 
the nonlinear propagator and examine if the memory of the cosmic fields to their initial conditions has the tendency 
to be kept by the non-Gaussianity or not. 

The organization of this paper is as follows. In section II, we review the basic equations of motion describing the 
evolution of the cosmic fields. In section III, we develop the general field theoretical framework so that the non- 
Gaussianity can be incorporated into the scheme. We have successfully calculated the generating functional which 
gives Feynman rules. We also develop the renormalization group method in this context. In section IV, we apply 
our formalism to calculate the nonlinear propagator for the non-Gaussian initial conditions and examine the effect of 
non-Gaussianity on the nonlinear scales. The final section is devoted to the conclusion. 



II. BASIC EQUATIONS FOR COSMIC FIELDS 

In this section, we review the standard Newtonian cosmological perturbation theory. Here, we consider the Einstein- 
de Sitter universe for simplicity. Of course, it is possible to extend our analysis to other cosmological models. 

First of all, let us consider the homogeneous cosmological background spacctimc. Taking the conformal time and 
assuming the flat space, we can write down the metric 

ds 2 = a 2 (r) [-fir 2 + ,\ ,!.r' . (2) 

The cosmological scale factor a is determined by solving FRW equations 

H = —^-a po , H = —a p , (3) 

where pg is the averaged density field and we have defined 7i = da/dr/a = a' /a. 

Now, let us consider the inhomogeneous distribution of the matter. The evolution of the total matter density is 
determined by the gravity including the effect of cosmic expansion. The actual density p(x, r) is deviated from the 
averaged density po(t) Let us define the density fluctuation as 

= p(x,r) - po(r) = f £ ik . x 
Po(t) J 

It obeys the equation of continuity and the peculiar velocity v is determined by the Euler equation in the presence of 
gravitational potential c/>. The gravitational potential itself is governed by the Poisson equation. Thus, equations of 
motion for the cosmic fields, S, v, and <f>, are given by 

g+ V -[(l + 5)v]=0, (5) 

— +Wv+(v- V)v = -V0, (6) 

OT 

V 2 0=^. (7) 

On large scales, the assumption that the peculiar velocity is irrotational would be valid. Then, defining 9 = V • v, 
we obtain the relation v(k, r) = ik#(k, r)/k 2 . After eliminating the gravitational potential, we obtain equations of 
motion in the Fourier space 



0(k, t)+ J d 3 qd 3 P 6 D (k q P)«(q, p)0(q, r)*(p, r) = 0, (8) 



dr 

^fc)+^(k,r) + ^(k,r) 



- J d 3 qd 3 p6 D (k - q - p)/3(q, p)0(q, r)0(p, r) = 0, (9) 
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where we have defined the Dirac delta function 3d and 



, * (p + q)-q a , \ (p + q) 2 p-q 
"(q, p) = 2 — ' Pii' p = — ^4 — ( 10 ) 



characterize the nonlinear gravitational coupling. 

It is now easy to solve the above equations perturbatively. Given the expansion 



*(k,T) = £* (n) (k>r), 0(k,r) = ][>W(k,T) , (11) 



the iterative solutions can be explicitly written down as functional of the initial fields Q. Once the probability 
functional for the initial fields are given, one can calculate expectation values of product of these fields. In fact, this 
standard formalism has been utilized to calculate the quasi-nonlinear evolution of the power spectrum, bispectrum, 
and other statistics such as genus statistics [H, H, 0j- In the nonlinear regime, however, further analytical tools are 
required to give more accurate results. This is crucial for the prediction of BAO feature, for example. In the next 
section, we introduce a useful approach to this end. 

III. PATH INTEGRAL FORMALISM AND RENORMALIZATION GROUP 

Now, we proceed to formulate the cosmological perturbation theory in the field theoretical manner. Doing so, we 
can utilize the idea invented in the field theory. In particular, the renormalization group method turns out to be 
useful. 

As is suggested by Crocce and Scoccimarro [8J|, the equations of motion can be rewritten in a convenient form by 
defining a two-component vector 

v f k r) = ( ^ l(k ' T) ^ - e-" ( 6{Kt) ) (12) 
M*, r) - y ^ (kj r) j-e y _ 9{K t)/h J , (12) 

where the index a = 1, 2 and 77 = loga(T)/a- ln denotes the e-folding number. The initial scale factor di n can be taken 
arbitrarily. We also define vertex 7 a {, c with 

7i2i(k,p,q) = i^(k + p + q)a(p,q) =7 U2 (k,q,p) , 7222^, p,q) = S D (k + p + q)/3(p, q) (13) 
and other components are zero. Using the above definitions, we obtain the equation 

(5 a bd v + ^afc) ¥>&(k, rj) = e ?) 7 afcc (k, -p, -q)p&(p, rj)<p c (q, rj), (14) 
where Q a b are components of the matrix 

" = ( -3/2 3/\ ) ■ < 15 ) 

Here, we have used the Einstein's sum rule 

<^ a (-k,0)^ b (k,0) = f d 3 fc(^ a (-k, 0)y> 6 (k,0) . (16) 

From now on, we should understand this convention is used when the same wavenumber vector appear twice in the 
same term. 

First, we consider the linear theory. The growing and the decaying modes can be written as u a and v a exp(— 577/2), 
respectively. Here, we have defined two basis vectors 

1 ) and Va = ( -3^2 ) ' (17) 
Let us define the linear propagator as 

(S a bd v + Qab) gbc(rja - r)b) = S ac S D (r/ a - rib) . (18) 
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Given the growing and decaying mode functions, we can write down the causal propagator g ab as 

gabiVa^Vb) = {u a Ub + v a Vbexp (-5/2(ri a -iib))}d{Va -Vb), (19) 
where u a and v a are dual vectors of u a and v a , 

u a = -{3,2) and fi = |(l,-l). (20) 
Note that we have the matrix components 

Uailb ^(n)' Vaib = \ ( -3 ~3 ) ' (21) 
Using this linear propagator, we can formally solve the equation (|14[) as 

Va(k,77 a ) = g ab (ri a ,0)(p b (k,0) + / drj b g ab (rj a , r] b )e n ^ bcd (k, -p, -q)<p e (p, JfoVdCq, ■ (22) 



Iteration gives the perturbative solutions. There are three building blocks, the initial field ip(k, 0), the linear propagator 
g ab , and the vertex 7 f, c - The initial field is usually assumed to have Gaussian statistics characterized by the linear 
power spectrum 

< ^(k,0>(k',0) >= (2Tr) 3 5 D (k + k')pt . (23) 

The nonlinear power spectrum can be calculated using the graphical method as was first demonstrated by Crocce 
and Scoccimarro Q. In particular, the renormalization method is utilized to perform the partial sum of diagrams. 
It turned out the approximation motivated by the renormalization method was quite successful [10j . In the next 
subsection, we consider more general statistics for the initial field (p a (k, 0) . 



A. Path Integral Representation 

What we are interested in are the statistical quantities characterizing the large scale structure in the universe. 
The statistics of primordial cosmic fields <p a (k, 0) are determined by the initial probability functional P[(p a (k, 0)]. To 
calculate desired quantities, we need to solve the nonlinear evolution equations and obtain the solution as a function 
of the initial fields. Namely, we have the statistics and the dynamics to be considered. More precisely, we want to 
calculate 

^expi J dr)J a (-k, r))ip a (k, r\\ <p(k, 0))\ 

= J dp a (k,0)Pfo> a (k,0)]expi J dr/J a (-k,r;)^ a (k,r;;^ a (k,0)), (24) 

where P[ip a (k, 0)] is the general probability functional for the initial field <p a (k, 0) and cp a (k, rj; ip a (k, 0)) is the solution 
of Eq. (fT4j) with the initial condition ip a (k, 0). This is a generating functional for multi-point correlation functions. 
Here, we shall combine the statistics and the dynamics in a unified framework. This can be achieved by the field 
theoretical path integral method. It is the path integral representation of the problem which can be used to perform 
the non-perturbative approximation. 

To derive the path integral representation for the cosmic fields starting from the expression (|24p . we introduce an 
auxiliary field tp a (k, 77) as 

J d<p a (k,0) J D<p a (k,r))P[ip a (k,0)]5 D (<p a (k,r)) - <^ Q (k, 77; (p(k, 0))) exp i J drj,J a (-k : ^(k, rj) , (25) 

where we used the Dirac delta function Sd- To separate the dynamics from the statistics, we use the operator L 
defined by 

L<^a = (Sabd v + fi afc ) (p b (k, r/) - e n j abc (k, -p, -q)<pi,(p, r?)y c (q, rj) ■ ( 26 ) 

Then, we get 



dip a (k,0) J Dtp a (k, rf)P[ip a (k, Q)]5 D (L<p a (k, rj) - tp a (k,Q)5 D (r))) expi J drj,J a (~k, rj)if a (k, 77), 



(27) 



where we have used the fact det L = 1 valid for the causal boundary conditions [19( . Here, it is convenient to introduce 
another auxiliary field Xa to exponentiate the Dirac delta function. This procedure is crucial to integrate out the 
initial field <p a (k, 0). The result becomes 

dtpaQs.,0) [ D<p a (k,r])Dx a {Kr))PlVa(k,0)] 



x cxp 



i / d?7X a (-k,77)L^ (k, 77) - ixa(-k, 0)Va(k, 0) +i / dr)J a (-k,rj)ip a (k,ri) 



= J Dtp a (k,Ti)Dxa(Kv) e 



C[xa(k,o)] 



exp 



' J d»7Xa(-k,77)L^ a (k,77) + i J drjj a (-k, rf)ip a (k, if) 



where we have defined the cumulant functional C[xo] by 

e c[x n (k,o)] = f ^ Q ( k5 0)P[^ Q (k,0)]e-^(-k,o)^(k,o)^ 



(28) 



(29) 



The field Xa(k, 0) is defined as the boundary value of the field Xa(k, 77). The cumulant functional completely charac- 
terize the initial statistics. It is possible to expand the cumulant as 

C[xa(k, 0)] = (2ir) 3 f-ixa(-k, 0)P a bXb(k, 0) + l -B abc (k x , k 2) k 3 ) x„(-ki, 0)x 6 (-k 2 , 0)x c (-k 3 , 0) + • 



(30) 

The first term is the power spectrum for the initial field. In the Gaussian case, only this term exists. In general, 
subsequent terms follow. The second term is called the bispectrum 

In the above action, the auxiliary field X a is introduced as a field. Hence, it is natural to add the source for this 
field. Thus, we have the final path integral expression for the cosmological fluctuations with the general statistics: 

Z[J ai K b ] = y > J D^ a (k,r ? )Dxa(k,r;)e c[x " (k ' 0)1 exp i J dry X a(-k, r^L^k, 7?) 

+i J d7]J a (-k,r])(p a (k 7 r]) + i J dr]K a (-k 7 r])xa(Kil) ■ (31) 

From the field theoretical point of view, C[x Q (k, 0)] can be regarded as the boundary action on the initial hypersurface 
rj = 0. In this sense, the field Xa(k, 0) is associated with the initial conditions. 



B. Generating functional 

Thanks to the auxiliary field Xa> we can perform the path integral completly. First of all, let us extract the nonlinear 
interaction part 

Z[J a ,K b ] = expiSintl-iJ-,-^} f D^ a (k,r 1 )D X a(Kv)e ClXa{K0)] 

O.J a 6K b J 



x cxp 



'J dr]x a {-k,r])Li iri (p a (k, r))+if drj,J a (~k, rj)(f a (k, r))+ij dr]K a (-k, rf)x a (K rj) 



, (32) 



where L; in is the linear part of L and Si nt denotes the interaction part of the action. As the field ip a is linear in the 
action, it is easy to integarte out it as 



6 S 

expiS lnt [-i — 7 -i—} I Dxa{k,i 1 )e clXa{k ' a)] 5 D {Xa'Lun + Ja)exp 
dJa oK b 



i / di]K a (k,r))xa(k,ri) 



(33) 



where XaX'Un = Xb(Vb)(~ d nb 5 ba + Q ba )- The constraint imposed by the Dirac delta function can be solved by using 
the equation 



as 



g a b(Va,Vb)(- d , lb S bc + fl bc ) = S ac S(r] a - r] b ) 
Xa(KVb) = - J drjJ a {k 1 Ti a )g ab {rj ai Ti b ) . 



(34) 
(35) 
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Now, it is straightforward to integrate out the field Xa- The final result is given by 

S S f 
Z[J a ,K b ] = expiS int [-i — ,-i-^-]exp C[- J dr]J a (k,r])g ab (r],0)] 



-i J dr)J a (-k,r]a,)gab('na,'nb)K b (k,r)b) . (36) 

This leads to Fcynman rules for calculating various quantities. Thus, we have shown that the graphical method is 
applicable to the evolution problem of cosmological fluctuations with the general non-Gaussian statistics. 

It is useful to see a more concrete expression for our path integral representation. In the case of the Gaussian 
statistics, the cumulant is determined solely by the power spectrum. Therefore, we have a simple expression 



Z[J a ,K b ] = cxp 



where 



dr]e v ^ abc (k, p, q) 
(2nf 



5K a (k) 



SJ b (p) 



<Wc(q) 



Zo[J a , K b 



(37) 



Z [J a ,K b ] = exp 



drj a dr) b J a (-k, rj a )P ab (r] a , rj b )J b (k, rj b ) 



-i \ dt] a dr] b J a (-k,r} a )g ab (r] a ,r] b )K b (k,r} b ) 



(38) 



This result coincides with the expression obtained by Matarrese and Pietroni. In the case of the non-Gaussian statis- 
tics, we have infinite series of the irredecible correlation functions. In this non-Gaussian cases, the extra contributions 
read 



Zo[J a , K b 



exp 



1 



( 2?r ) ~o / dr) a dr) b J a {-k,rj a )P ab (r] a ,r] b )J b (k,T] b ) 



-g / dr} a dr} b dri c Ja{-k,r}a)Jb{-K'n b )J c {-k,r) c )B abc (ki,k 2 ,k 3 j 



-i J di] a dr/ b J a (-k,r) a )g ab (r] a ,r) b )K b (k,r] b ) , (39) 
where we have defined the propagated initial correlation functions 

P ab(Va,Vb) = gac(Va,0)g b d(Vb, 0)P cd , (40) 

B a bc (k 1 ,k 2 ,k 3 ,T] a ,T] b ,ri c ) =g ad (r) a ,0)g be (r) b) 0)g c f(r) c ,Q)B de f(ki,k2,k 3 ) . (41) 

The formula ()37|) gives Feynman rules. In this paper, we will be concentrated on the effect of the initial bispectrum 
B abc - 

In the quantum field theory, the generator of the connected Green functions can be written as 

W[J a ,K b ] - -i\ogZ{J a ,K b } . 
The expectation values 4> a and Xa can be defined as 

SW SW 



Motivated from the linear expression 



(27r) 3 P L 6 = -i 



we define 



(27r) 3 P afc = -i 



S 2 W 
5J a SJ b 

S 2 W 

SJ a SJ b 



OJa 0K a 

S 2 W 



(42) 
(43) 



,9ab — 



■, Gab — ~ 



5J a SK b 
5J a SK b 



(27T) 3 Sf fcc 



, (2ir) 3 B abc 



S 3 W 
6J a 5J b SJ c 

5 3 W 
SJ a SJ b 6J c 



Carrying out the Legendre transformation of W, we can define the effective action 

r[^a,Xa] = W[J a [ip a ,Xa],K a [Lp a ,Xa}} ~ J dT]d 3 k(J a [ip a , Xa}fa + K b [(fi a , Xa]Xb) ■ 

where 

Ja[<Pa,Xa\= —TT , J<a[fa,Xa\ = -J • 

0(f>a OXa 



(44) 
(45) 
(46) 
(47) 



a 



9ab 



+ 



- zexp(?7)7 Q 6 



a 



a 



pL 

^ ab 



tdL 

D abc 



FIG. 1: Feymann diagrams; an arrow represents the direction of the time. 



C. Renormalization Group 



The idea of the renormalization group introduced by Matarrese and Pietroni is as follows. First, we introduce a 
filter function with the UV cutoff A in the P° b and B® bc . This defines a fictious theory where the linear perturbation 
theory works well. We denote various quantities in the cutoff theory with the suffix A, for example like as P a b,x- 
When this cutoff scale A goes to infinity, the original theory is recovered. As A becomes large, nonlinear effects arc 
incorporated gradually. This process can be expressed by the renormalization group equation. 

The RG equation for the effective action can be deduced as follows. Taking the derivative of the generating 
functional Z\ with respect to A, we obtain 



d x Z x = 



(2tt) e 



drjadr] b d\P ab ,\(p)6(r] a )6(ri b ) 



s 2 z x 



5K a (-p,r])SK b (p,r)) 



(2nf 



dr) a dr) b r] c d\B abCt x(pi,P2,P3)fi('na)5(r)b)S(r) c ) 



s 3 z x 



SK a (pi,r])SK b (p2,r])SK c (p3, rj) 



(48) 



where we have replaced the auxiliary fields \a by the functional derivative with repect to K a . The above equation 
can be translated to the equation for the generating functional of the connected Green functions. Using the relation 
(|42l). we obtain 



d x W x = 



(2^ 



drj a dr] b d\P ab ^5(r] a )6(r) b ) 



S 2 W X 



+ i- 



SW 



sw 



6 



8K a (-p, rj)5K b (p, rf) SK a {-p, rf) SK b (p, rf) 
dr] a dr] b dr] c d x B abC:X (p 1 ,p 2 ,p 3 )S(r] a )S(ri b )S(r] c ) 
S 3 W ( SW 



S 2 W 



SK a {p ll r))5K b {p 2l r))5K c {p 3l r)) \5K a (pi,r)) SK b {p 2 , rf)5K c (p 3 ,rf) 
SW SW SW 



+ (eye. [a, b, c}) 



SK a {p ll rf) SK b (p 2l r)) SK c (p 3 ,r l ) / 
Now, we can write down the RG equation for the effective action. Using the definition (|46|) . we get 

= d\W\[J\,a,[tp a ,Xb],K\,a[Pa,Xb]] ~ (dx Jx,a[fa, Xa]fa + d X K X , b [if a , Xa]Xb) 



(49) 



(2n) s 



di] a dr] b d\Pa b ,\S(r] a )S(ri b ) 



5 2 W X 



+ t- 



SW 



SW 



6 



8K a (-p, rf)5K b (p, rj) SK a (-p 7 rf) SK b (p, rf) 
dr) a dr] b dr] c d\B abC: x(pi,P2,P3)S(r)a)S(r] b )S(ri c ) 

S 3 W . ( SW S 2 W 



+ 1 



8K a {p ll rf)5K b (p 2 ,rf)SK c (p 3l rf) \5K a (px,r)) 8K b (p 2l rf)SK c {p 3l rf) 
SW SW SW 



+ (eye. [a, b, c]) 



^a(Pi,7?) SK b {p 2l rf) 5K c (p 3l rf) 



SW SW 

+ -fl—d\J\,a + : Ti r d\K X ,a ~ (d X J X . a [f a ,Xa}fa + d X K X , b [<Pa, Xa\Xb)- 

oJ\ a SK\ a 



(50) 
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The last line vanishes due to the definition of the expectation values. It is convenient to separate out the contribution 
of the tree part from the effective action. Let us define the tree part 



To : a = / drj a di] b 



i(2nf 



Xa(-P, v)Pab,x{q)S{r]a)S(T]b)xb(p, rj) + Xa{-p, ri)g a b<Pb{p, i) 



+ 



(2tt) S 



di] a dr] b dr] c B abcA (p 1 ,p2,P3)S{r]a)S(r] b )S(T] c )xa(~Pi,v)Xb(~P2, v)Xc(-Pa, v) (51) 



and the interaction part 



int.X 



0,A 



(52) 



Then, Eg. ([50)) can be written as 



d\Tint,\[<Pa,Xb] 
(2nf 



dil a dr]bd\Pab,\d{Va)fi(Vb) 



5 2 W X 



5K a (-p,r))6K b (p,r)) 



(2tt) s 



df] a df] b dr] c dxB abCt x(pi,P2,P3)5(Va)5{Vb)HVc) 



S 3 W 



SW 



S 2 W 



8K a (pi,rj)8K b (p2,rj)8K c (p a ,rj) 8K a (pwij) 8K b (p 2 , r))SK c (p 3 ,r)) 



+ (a <-> b «-> c) 



(53) 



Ja = -f-,K a = JL- 



This is the RG equation for general Green functions. In the next section, we will apply our formalism to the calculation 
of the nonlinear propagator. 



IV. NONLINEAR PROPAGATOR IN THE PRESENCE OF NON-GAUSSIANITY 



Nonlinear interaction causes a deviation from the linear propagator g ab . This effect is interpreted as a propagator 
renormalization by Crocce and Scoccimarro [18| . Physicaly, it is interpreted as a measure of the memory of initial 
conditions. The nonlinear propagator is defined by 



G a b(k,T]) 



< p o (k,0)<p 6 (-k,0) > 



(54) 



If the linear approximation is good, the above expression should give the linear propagator g a b(r])- 

We shall calculate the nonlinear propagator by employing the renormalization group equation. Doubly differenti- 
ating Eq. (|4l))) with respect to K a and J a , we can get the RG equation for the propagator as 



s 2 w x 



8J a (-k, r] a )8K b (k', rj a ) 



-8(k - k')d\G ab ^(k] -q a r] b ) 



di] c d7] d dxP cd}X S (r) c )8 (i] d ) - 



S 4 W, 



8J a (-k, r) a )SK b (k', r] a )SK c (-p, r] a )SK d (p, r} a ) 

dr] c dr] d dr] e d\B cdet x(pi,P2,P3)8{r] c )S(T] d )8(r] e ) 

S 5 Wx 

SJ a (-k,r] a )SK b (k' ,'q a )SK c (-p 1 ,^ a )SK d (-p 2 ,r]a)SK e (-p 3 ,r] a ) 



(55) 



d/d\—$- 




FIG. 2: RG equation diagram for the propagator. 



In the above, we have used the fact 1 

S n W[J a ,K b ] 



„ = , for any n. 



5K ai {kum) ■ ■ ■ SK an (k n ,ri n y J ' =0 ' K ' 
Now, we shall rewrite connected correlation functions by the irreducible correlation functions 

,(n) _ 8 n T[ip a ,Xb] I 



5 D (ki 



S Xai (h,rn) ■ ■ ■ 6<p an (fc„, rin) 1 va=o.xo=o - 

For example, the four-point function can be written as 

S 4 W X 
8J a 6K b SK c SK d 

dsi ■ ■ • ds^G aej x{ki;r] a , si)Gfb,x(k 2 ; s 2 , ri b )G gc .x(k 3 ] s 3 , r} c )Ghd,x(^ s±, r/d) 
:{- J ds 5 ds 6 Gij,\(ki - k 2 ;s 5 ,s 6 ) 

xT xIv>fV*( kl > Sl > fc2 ' S2 > kl ~ fc 2,S5)rg ) ¥ , jVh (fci - k 2 ,SQ]k 3 ,s 3 ;k 4 ,Si) 

+ r XeWv s Vh( fcl ' Sl ! fc 2,S2;fc3,S3;fc4,S4)}^(fcl ~ k 2 ~ k 3 - fc 4 ). 



(56) 



(57) 



(58) 



In Figl^l we have drawn the Feynman diagrams for the RG equation. The four-point function in the RG equation 
gives the first and second diagrams in Figj2] Other diagrams in Figj2] come from the five-point function 

6 5 W X 



SJ a 5K b 6K c 5K d 5K e 

dsi ■ ■ ■ dsiG a f } x(ki;T] a , si)G gb ,x{k 2 ; S2,r)b)Ghc,\{kr, S3,r]c)Gict,x(k4] s<t, r)d)Gj e ,x(k5; S5,r] e ) 



x{ I ds 6 ds 7 ds s ds 9 GkLx(ki - k 2 ; s 6 , s 7 )G mn ,x(h + k±] s 8 , s 9 )T < £^^ k (k 1 , si; k 2 , s 2 ; k x - k 2 ,s 6 ) 

xT^ VhVm (ki - k 2 , s 7 ; k s ,s 3 ; h + fc 4 , s s )T^ VzVj (k 5 + fe 4 , s 8 ; fc 4 ,s 4 ; k 5 ,s 5 ) 

+ ...}S(k 1 -k 2 -k 3 -k i -h) , (59) 



1 This can be proved as follows. First of all, we need to keep it in our mind that the propagator is causal. Suppose that we take 
derivatives of W with respect to only K a (not J a ), then propagators of all external lines go to the future direction. Now, let us consider 
the propagator from the most future external line. This propagator must connect to a vertex. However, one of the edges of the vertex 
still directs the future. This means that the external line is needed in the future, which contradicts the initial assumption. 
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where we have explicitly displayed only the relevant one. 

We can now calculate the RG equation. On the nonlinear scale, modes enter into the damping phase. Therefore, 
the nonlinear effect allows us to consider only the low wavenumber mode for the internal lines. Hence, we can take the 
large k limit, where k is the wavenumber of the rcnomalizcd propagator. In this paper, we approximate a solution by 
taking into account only the running of the 2-point function. We also ignore the propagation of the decaying mode. 
Thus, the coupling constant ^ a bc always appears in the following form: 



7(k,P,q) 



Ualabc{K P, q)u b U c 

IcW (6pV + 5(p2 



q 2 )p-q + 4(p-q) 



In the limit p 3> q, we obtain 



7(k,P,q) 



p • q 



(60) 



(61) 



In our approximation, the diagrams we need to consider are the first, third, fourth and fifth ones in FigEJ It is known 
that the fourth and fifth ones are higher order contributions than the third one [18| . The third one directly interacts 
with the initial conditions, so that all informations of the intial conditions transmit. From Eq. (|60[) . we see the direct 
interaction linearly grows with the wavenumber p. On the contrary, there are interactions outside the main path in 
the fourth and fifth diagrams. There, we can not take the limit k ^> q and, hence, we cannot expect the linear growth. 
Thus, we can ignore them. As an additional approximation, we replace the propagator in the right hand side of the 
RG equation in Fig[2]with its linear expression. 

Using the above aproximations, we can get the RG eqation for the propagator, 

d\G abi x(k;rj ai r] b ) 

= 9ab(r]a,rib)\ 2 2 (27r) 3 / ds 2 / dsi exp(si + s 2 ) 



x J d 3 qd x Pcd,x(l) il cUd-fOi 1 q, k - qh(k - q, -q, k) 



Va. 



S3 r s 2 



-2 3 (2?r) 3 / ds 3 / ds 2 / dsiexp(si +s 2 + s 3 ) 

J rib J rib J r/b 

x J d 3 q 1 d 3 q 2 d 3 q 3 dxB cde .x(q_ 1 ,c[2, q.z)u c u d u e 

x 7( k , qi,k-q 1 )7(k-q 1 ,k + q3,q 2 )7(k + q3,q 3 ,k) 

= gabiVa, Vb)^2 (exp(?7a) - exp^f,)) 2 d\F Pt x[k] - 8 (exp(?7 Q ) - exp(?7 6 )) 3 d\F B .\[k 
where we have defined 

e (2tt) 3 J d 3 qP cdjA (g)u c u d 7(k,q,k-q)7(k-q, -q, k) 



(62) 



(2n) 3 ^ f dqP cd Aq)u c u d (^ ~ ^ + 9, 2 - \* 



3fc 5 9fc 3 9 21 q 3 9q 5 . 



k + q 



k — q 



(63) 



and 



F B ,\[*] = 



(2?r) 3 J d 3 q 1 d 3 q 2 d 3 q 3 B cde ^((li,(l2,(l3)u c u d u e 

X7(k, qi, k - qi)7(k - qi, k + q 3 , q 2 )7(k + q3, q3, k). 



(64) 



Here, we promote the linear propagator to the full nonlinear one following Mataresse and Pietroni [12j. Solving the 
differential equation (f6"2"|) with the initial condition G ab .\ = o — g a b and taking the limit A goes to infinity, we obtain 



G a b{k; ri a , rj b ) = g a b{k; r} a , rj b ) 

x exp|2 (exp(r/ Q ) - exp(?7 6 )) 2 F P [k] - 8 (exp(r? a ) - exp(ry b )) 3 F B [k}\ 



(65) 
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As a concrete example, we consider the model 

^o(x) = % (x) + W^(x) - <^}), (66) 

where t/? g (x) is a Gaussian held. It should be noted that this £nl is different from Jnl, because we are considering 
ip a instead of the curvature perturbation <F Then the initial bispectrum B a \, c is written as 

B abc {q ai q b , q c ) - 2^ NL S(q a + q b + q c ){P(q a )P(q b ) + P{q b )P{q c ) + P(q c )P(q a )), (67) 

where P(q) = P c d(q)u c Ud- In the k^> p limit, we can deduce 

F P [k] ► -(2tt) 3 ^ f dqP cd (q) , (68) 

and 

F B [k] ► (2tt) 3 ^- f d 3 qi d 3 q 2 d 3 q 3 

^ (k-qxKk-qa) (k.q 2 )(k.q 3 ) (k ■ q 3 ) (k • qQ \ 

x 31 + 31 + 312 £(qi,q 2 ,q 3 ) , (69) 

V 12% IsHl / 

where we have defined B(q) = B c d e (q)u c UdUe- 

From the shape of the CDM spectrum, we can estimate the scale where the running of the propagator is important. 
The nonlinear effect works when the argument of exponential in Eq. ((65|) becomes 0(1). The threshold for the power 
spectrum (Fpclq) becomes k ~ 0.1(1 + z)Mpc~ . On the other hand, the threshold for the bispectrum (FgOg) is given 

by k ~ 0.03£~|(1 + z) 3 / 2 Mpc _1 . This means that, if the statistics is positively skewed and £jvl > 0(1), the effect of 
the initial bispectrum is a dominant one at present in the non-linear propagator. On the other hand, the negatively 
skewed non-Gaussianity delays the time each comoving mode enters the nonlinear regime. The above arguments 
imply that the non-Gaussianiy significantly affects the running of the propagator. 

The propagator can be regarded as the measure of the memory of initial conditions. Only in the case the memory is 
left in the data, we can detect the primordial non-Gaussianity. According to the behavior of the nonlinear propagator, 

the modes in the range k > 0.1(1 + z)Mpc 1 or k > 0.03^(1 + z) 3/2 Mpc 1 seems to have already lost the memory 
to the initial conditions. This means that, at the time z ~ 1, we can not observe the primordial non-Gaussianity 
around the scale k — lMpc^ 1 . In order to investigate the primordial non-Gausianity at the comoving scale IMpc, we 
have to look at the Universe at z > 10. One possibility is the observation of 21cm line. As to this, many observational 
projects, such as the Low Frequency Array (LOFAR), the Square Kilometer Array (SKA), the Mileura Widefield 
Array (MWA) and the 21cm array (21CMA), are ongoing or being planed [3l| . 

V. CONCLUSION 

We have studied the Newtonian cosmological perturbation theory from the field theoretical point of view. We 
have extended the renomalization group method proposed by Mataresse and Pietroni to the case of a primordial 
non-Gaussian density fluctuations. In particular, we have obtained the generating functional for the cosmic evolution 
of fluctuations with non-Gaussian statistics. As an application, we have calculated the nonlinear propagator and 
examined how the non-Gaussianity affects the memory of cosmic fields to their initial conditions. It turned out that 
the initial non-Gaussianity affects on the running of the propagator. For the positively skewed case, the nonlinearity 
starts at ealier stage. In the opposite case, the nonlinearity is postpond compared with the Gaussian case. 

Assuming the 0(1) positively skewed non-Gaussianity, we can conclude that the nonlinear propagator damps in 

the range k > 0.1(1 + z)Mpc _1 or k > 0.03^^(1 + z) 3 / 2 Mpc~ 1 . Hence, if we want to investigate the initial non- 
gaussianity at the scale k > IMpc -1 , the observation at z > 10 is required. One interesting possibility is to observe 
fluctuations of 21cm absorption line in the cosmic maicrowave background radiation. In the negatively skewed case, 
the effect is opposite. Therefore, the time entering the nonlinear scale is postponed for the mode with a fixed comoving 
wavenumber. 

As an application of our formalism, we can estimate the effect of the primordial non-Gaussianity on the BAO. More 
intriguingly, we can calculate the bispectrum of cosmic fields which provide a more clear test of non-Gaussianity. It 
is also intriguing to calculate BAO feature in the bispectrum. Again, the bispectrum of 21 cm line fluctuations is an 
interesting target. 

Recent observational progress allows us to know the large scale structure at high redshift. It should be emphasized 
that the field theoretical approach gives a simple way to calculate correlation functions with different times. Therefore, 
the field theoretical method in cosmology deserves further investigations. 
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